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Abstract 

We show how to compute real space renormahzation group flows in lat- 
tice field theory by a self-consistent method which is designed to preserve 
the basic stability properties of a Boltzmann factor. Particular atten- 
tion is payed to controlling the errors which come from truncating the 
action to a manageable form. In each step, the integration over the fluc- 
tuation field (high frequency components of the field) is performed by a 
saddle point method. The saddle point depends on the block-spin. Higher 
powers of derivatives of the field are neglected in the actions, but no poly- 
nomial approximation in the field is made. The fiow preserves a simple 
parameterization of the action. In the first part the method is described 
and numerical results are presented. In the second part we discuss an 
improvement of the method where the saddle point approximation is pre- 
ceded by self-consistent normal ordering, i.e. solution of a gap equation. 
In the third part we describe a general procedure to obtain higher order 
corrections with the help of Schwinger Dyson equations. 

In this paper we treat scalar field theories as an example. The basic 
limitations of the method are also discussed. They come from a possible 
breakdown of stability which may occur when a composite block-spin or 
block variables for domain walls would be needed. 
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1 Renormalization group flow on the lattice 



Suppose we start from the Euclidean action S'°[(/i''] of a field theory which lives 
on some lattice A'^ or on the continuum. It depends on a field (fp on hP . We 
wish to compute a sequence of effective actions 5* which live on lattices A* 
of increasing lattice spacing a^, i=0,I,2,... . 

In principle the sequence of actions is defined once we fix a block-spin def- 
inition. The block-spin definition determines a field $ = 0*"'"^ on the coarser 
lattice A'+^ as some kind of average 

$ = cv 

of the field (p = (p^ on the lattice A*. The actions are defined recursively by 

= e-^"^'[*l = J D(t)S{C(t) - $)e-^'['^l. (1) 

For given is the partition function of an auxiliary theory in which 

only those variables are integrated out which are to be interpreted later as high 
frequency modes of the field. 

We will also need certain expectation values of this auxiliary theory. If 
A ^ A[(p] is an observable, we set 

{Ay = J L>0A[0](5(C> - $)e-^'['^l. (2) 

These expectation values depend on but we neglect to indicate this depen- 
dence explicitly. In the calculations, two auxiliary <i>-dependent quantities will 
play an important role, the background field ^ '[$] and the fluctuation propagator 
r*[vli*]. The fluctuation propagator is considered to depend on the block-spin $ 
through the background field — 5"* [<&] . Both auxiliary quantities are defined 
as expectation values of the auxiliary theory. 

W^{z) EE ^^[<i>]{z) ^ Wz))\ (3) 

rmiz,w) = wz)^{w)r-{^{z)r{^{w)y. (4) 

These definitions differ from those of Gawedzki and Kupiainen, cp. In 
their work, expectation values of a free theory were used. The present approach 
is more in the spirit of Balaban's work on gauge theories, where expansions of 
the full action around its minimum were used for reasons of gauge invariance 

i- . . 

The validity of the method requires that the auxiliary theory has an infrared 
cutoff of order a~[_\ . Therefore the fluctuation propagator must decay exponen- 
tially with decay length no larger than about a^+i. When this condition is 
violated, it signals a bad choice of the block-spin definition. 

We are not able to calculate the functional integrals (|l],^ exactly. There- 
fore approximations are necessary. The calculation is done in a self-consistent 
way. Both the background field and the fluctuation propagator are used in the 
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calculation. Conditions will be imposed on them which ensure that equations 
(1) and (|) arc fulfilled within the accuracy of the approximation. We discuss 
in subsection |2.4| how a truncation of the effective action to a manageable form 
can be justified. No expansion in powers of fields is involved. 

It follows from the definition (^) that the fluctuation propagator satisfies the 
constraints 

Cr['^']^0 = r[^']C". (5) 

The simplest choice of block-spin for a theory of scalar fields is as follows. 
One identifies the sites x of the coarser lattice A*+^ with disjoint hypercubes in 
the lattice A% and one chooses the block-spin as the average of the original field 
on these hypercubes 



^{x) = / Cix,z)(j){z) = av0(z). 

Jz 

One can try to improve on the locality properties of the effective actions by 
modifying the block-spin procedure p^ . 

We may imagine starting from a lattice of finite volume. After a finite 
number of steps we arrive at a lattice A^ which consists of only a single site. 
The field <i> on this lattice is some average of an average ... of an average of the 
original field 0°. Let us interpret it as magnetization. The action 5^[<i>] will 
give us the constraint effective potential - i.e. the free energy - as a function of 
the magnetization. 

We will make the following assumptions on the choice of block-spin in this 
paper. We assume that is a linear map, and that C*C" is invertible. (For 
gauge theories, linearity will have to be given up.) We assume also that C" 
interpolates constant fields to constant fields, and we impose the following nor- 
malization condition 



C>(x) = / <^(z) 

for constant fields. For the above mentioned choice of block-spin, the equation 
is an identity for arbitrary fields 0. 

Let W be the Hilbert space of square summable functions on A*, and let 
be the subspace of functions C with vanishing block average C*C' — 0. 
Invertibility of C'C" ensures the existence of the orthogonal decomposition 

where C = {rC")~'^C\ This is true because every (jj € Ti.^ can be decomposed 

as ^ 

= C"$ + c 

with $ = e and = </> - C"<i>. The first summand is obviously 

in C"7i'+^ and the second is in Hq because it vanishes when we apply C: 
= C'0 - C*C"<i> = $ - $ = 0. Also the orthogonal decomposition is true 
because the scalar product vanishes: (C,C"$) = {{C'C"y^CX\ *) = 0. 

Notation: We will use letters z, w... for sites in A* and x, y... for sites in A'+^. 
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2 The saddle point approximation 



2.1 Parameterization of the action 

We wish to compute the action S'*+^[<i>] which depends on a field $ on the 
lattice A*+^ from the action S'^[(j)] on the finer lattice A*. We will make some 
approximations to perform the calculation. We present a parameterization of 
the actions which will be preserved by the approximate renormalization group 
flow. 

The action S'*+^[<i>] will depend on the field $ through a functional 
The field ^E"* lives on the lattice A*. It is called the background field. 
The recursion formula will have the form 

= - itrlnr[**], (6) 

^ji = ^^[$]. (7) 

Approximations will be made such that 

1. The parameterization of the action preserves its form 

2. The fiuctuation propagator, the effective action and its first two field 
derivatives, and the background field ^''[$] and its derivative \E'*'[<i>] with 
respect to the field $ will enter into the recursion relations, but they need 
only be calculated for constant fields. This computational problem is fit 
for a PC. 

We will use a dual notation for derivatives with respect to fields 

and similarly for second derivatives ^zw 

The background field 5** will be determined as a functional of $ by the 
saddle point condition which involves the previous action 

= min subject to C*** = (f> . 

The fluctuation field propagator F* lives on lattice A\ It is a selfadjoint 
integral operator with kernel 

r'[0](z, w) , z,w e 

and is defined as pseudoinverse of S** "[</>] for all </> in the sense that 

r[<j)]s"'[<j)]r[(i)] = r[<j)]. (8) 

In the case where the field (p is equal to the background field [$] , the fluctu- 
ation field propagator is positive definite on W^. It satisfies the constraints (||), 
viz. 

c'r[-^'] = o = r[^']c". (9) 
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Therefore the space W is mapped into by F*, and vanishes on C"?i*"'"^. 
The trace tr which appears in the action (^) is to be understood as a trace over 
tic- 

We will explain later on how the background field and the fluctuation propa- 
gator can be evaluated by a self consistent approximation which involves neglect 
of higher order terms in gradients of fields. 

It is the crux of any real space renormalization group method to find trun- 
cated forms of the effective actions which can be parameterized in a manageable 
form. 

Often this is done in an ad hoc fashion which throws away some pieces which 
are irrelevant in a perturbative sense - higher powers of fields, for instance. This 
is not really justified because existing irrelevant terms get suppressed in the next 
RG-step. But new irrelevant terms are created by the marginal and relevant ones 
at the same time. As a consequence, there is a kind of equilibrium, so that (along 
the renormalized trajectory) the irrelevant terms are determined by the relevant 
and marginal ones. They are not necessarily very small and they influence the 
flow of the marginal and relevant coupling constants. Throwing them away after 
each renormalization group step introduces therefore a systematic error which 
accumulates. 

In principle there is a better way. One could determine the irrelevant pieces 
as a function of the marginal and relevant ones by solving a fixed point equation. 
But until now this is practical only in simplified (hierarchical) models . 

The motivation for our method of truncation will be described below. We will 
argue that our approximation becomes more and more accurate the larger the 
scaling factor (ratio of lattice spacings) s ~ at+i/ai. We will present numerical 
evidence which confirms this. 

On the other hand, the saddle point approximation becomes exact in the 
limit s I— !■ 1 (when the nature of the cutoff permits such a limit). It becomes 
less accurate with increasing scaling factor because the phase space for high 
frequency modes increases. 

The truncation consists of a local approximation to the second derivative 
5'" = —W oi the action. For more precise notation, set 

-S\,4<l>] - W'mz,w) (10) 
{W'[q^]f){z) = / Wmz,w)f{w) (11) 

We will approximate by a function iy(z, w|^) which depends only on the 
value ^ of the field at one site of the lattice A* . If / G W is a function on A* 
then we approximate 

{W'W){z)^\j \^\zM4>{z)) + W'{zM4>{w))\f{w) . (12) 

can be determined from the knowledge of for constant fields 0. Sym- 
metrization is performed to maintain hermiticity of W . 
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To understand the meaning of the approximation, suppose that is the 
action of a 0^-theory. Then 



{S^"f) (z) = 



/W- (13) 



We see that our approximation is exact in this case because the field dependent 
part of S**" only involves the field at a single point z. 

In particular, the approximation is exact for free field theories. Since the 
saddle point approximation is also exact for free fields, our method is exact for 
free field theory. 

A general action can be decomposed into a potential and a term whose first 
derivative with respect to the field vanishes for constant fields. We call it a 
generalized kinetic term. 

= generalized kinetic term + V''[(f)] . (14) 

A field independent contribution to is of no interest. Therefore it suffices to 
know the derivative V^' . It is uniquely determined if we know the derivative S^' 
of the action for constant fields. 

Our method consists in deriving recursion relations for V'^ and i.e. for 
the first and second derivatives S"' and W"^ = — S**" of the action evaluated at 
constant fields. 

Moreover we wish to obtain the effective actions as a function of the block- 
spin and not only as a function of the background field. Therefore we will need 
information on the block-spin dependence of the background field. It turns out 
that the use of the recursion formulae requires knowledge of the derivative of 
the background field with respect to the block-spin. There is a formula for this: 

^\Am^) ^ ^^*''[*] = (i-r'snc^*{z,x). (15) 

For a block-spin $ which lives on a lattice made of a single site, translation 
invariance implies that there exists a saddle point ^'^ = 'I'*[<I>] of the action which 
is also constant. Because of the constraint C'^' — $ and our normalization 
conventions it follows that 

^f'[$]^C"$. (16) 

So we know the background field in this case, assuming the saddle point is the 
minimum. If it is not, this signals a breakdown of stability. We will come 
back to such possibilities. They are familiar from the Maxwell construction in 
thermodynamics . 

The argument can be extended to show that to any block-spin $ which is 
constant on A*"*"^ there exists a saddle point (|l6|) which is constant, assuming 
that translational invariance is not spontaneously broken. If the invariance 
under translations by block lattice vectors is not spontaneously broken in the 
auxiliary theory with expectation values (0) , then the appropriate saddle point 



6 



will be invariant under block lattice translations. Therefore it is equal to 
the constraint minimum of the action on a single block with periodic boundary 
conditions. Now we can appeal to translation invariance of the action under 
translations on the fine lattice A* again to conclude that there exists a constant 
saddle point unless translational symmetry is spontaneously broken. It may 
happen that it is spontaneously broken or the saddle point is not a minimum. 
Again, this signals a breakdown of stability. 

In view of its interpretation as a propagator for the high frequency modes 
of the field which are integrated out in a renormalization group step, the kernel 
r'[(/)](z,w) of the fluctuation field propagator must decay exponentially with 
distance \z — w\ with decay length at most one block lattice spacing Oi+i (see 
. This is a condition which limits the range of applicability of the present 
method. It could be monitored during the computation of the renormalization 
group flow. When it is violated, this will typically entail a violation of the 
locality properties of the next action 5*+-^ as well. Locality requires that ^zw 
decays exponentially with the distance between z and w with decay length no 
larger than one lattice spacing ai. 

When this is violated it is typically a sign that the choice of the block-spin 
definition is not appropriate. This can happen for several reasons, see section 

I 

2.2 The detailed procedure of saddle point approximation 

Let S*' be a generic action. We now want to evaluate the effective action 5*+^ 

g-s'+i[<i>] ^ j r)^s{C4> - $)e-^'['^] (17) 

in a saddle point approximation. For this purpose we split the field 

(/)(z) = *^[$](z) + C(z) (zeA^). 

into the background field '^^ [<&] and a fluctuation field C ■ The background field 
has also the meaning of a mean field, because in the Gaussian approximation 
that we will use for the action eq.(^ holds 

The mean field 5"* is determined as a functional of the block-spin field $ by 
a saddle point condition which involves the previous action 

S'*[vl''] = min subject to C**' = (f> . 

The result is a nonlinear equation for [$] 

S^'l^S'] = C"X\ (18) 
C'*' = $ . (19) 

A'[<I'](a;) are Lagrange multipliers, x £ A*+^. 
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Expanding the action around the mean field to second order in C yields 

Jz ^ Jz,w (-20) 

The linear term in the expansion vanishes 

{s",C) = (C"A^C) = {x\cc) = 0, 

smce due to (|l|) CX' = 0. Therefore we get 

g-s-+i[*] _ 1 £,^'-^(c(»)e--5'[*'[*l]-iC''5"'[*'[*]]C' (21) 



n 

s'[*'[i>]] / n/^v-K'r'[*'[i>]]"'C 

H 



^^,^_iC'r-[*'[*]]-^C' (22) 



Note that S"^" is a linear operator on 7i*. In appendix A.l we show that can 
be written as 



V' = S"' ^-S'" ^C"{CS"' ^C"y^CS'" \ (23) 



Because of CT* = ~ PT** the fluctuation propagator vanishes on C^*7^ 



and is a linear operator W TCq. Therefore we can use the quantity (F*) ^ 



o/,; in (22). The extension of F' to all of W is a pseudoinverse of 5"' as 
described by (|8|). 

The remaining fluctuation integral is then Gaussian with covariance F* [5** [$]] 
and may easily be performed. The result is 

S'»+i[$] = S'*[5f*[$]] - itr In F^ (24) 

The trace tr is to be understood as a trace over . 
If one uses in (21) 

6{4>) ^ lim e"*/."^^'^' , 

K — >00 

one obtains the following alternative expression for the fluctuation propagator. 

F*[1'*[$]]= lim (5*"[**[$]] + kC"C^)"\ (25) 

K — >00 

2.3 Recursion relations 

By iterating the saddle point approximation we obtain 5*^^ from S^^^: 

S-'+^M = S'*+i[^'''+i[(^]] - itrlnF'+i[1'''+i[(^]]. 

Now we need F*+^. It is defined in terms of 5"*+^" by an equation similar to 
(p3|). Therefore we try to get a recursion relation for 
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In appendix |A.5| the following recursion relation will be deduced by differenti- 
ating eq.(|^). (We neglect to write functional '5' [$] -dependencies on the right 
hand side.) 



z,w£A' 
1 

2 



r{zi,Z2)W\4z2, Z3)r{z3, Zi)W\^{zi, zi)*\,(z)^\j,(«;) 



+ ^ J r{zi,Z2)W\Uz2,z,)^\,{z)^\y{w) (26) 



Z.W,Zi .Z2 GA"^ 

[ 

2 



+ i / r\zi,Z2)W\,{z2,Zi)riz,Z3)W\Uz3,Zi)^\yiw)^\.Az4)- 



This functional relation is an exact consequence of eq.(^, no localization 
approximation has been made yet. 

We see that we need an equation for the derivative ^f*' of the background 
field. In appendix A. 2 it will be shown that 

= (1 -r5'")(?*. (27) 

The recursion relation is shown in graphical form in figure ^. The graphical 
notation is described in figure |l|. 

We see that the recursion relation for S**" contains a 1-particle reducible 
term. This is an artifact of the simple saddle point approximation which will 
disappear in the improved version of section ^ 

As described in section 2.1 the action 5* can be split into a potential term 
and a generalized kinetic term whose first field derivative vanishes for 
constant fields; the potential term is determined by the derivative S"^' evaluated 
at constant field. 

We give the recursion relation for the first derivative 5*' in a general form. 
Using (see appendix A. 3) 

Y\^=V^W\X ■ (28) 
one gets by differentiation of ( p^ 

^ ^ (29) 

The graphical illustration is shown in figure 0. 
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S(l){zi) ■ ■ ■ 6(j){Zn) 



Zi - ■■ z„ 



Q 

• • • 



(5$(x) 



Figure 1: Definition of graphical notation for the derivatives of the action and 
of the background field 



^^^^^ = <wwv«^^B>/ww» + 




Figure 2: Recursion relation for = — 5'" in saddle point approximation 



Figure 3: Recursion relation for —5*+^' in saddle point approximation 
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For constant block-spin field the derivative of the generalized kinetic term 
vanishes. Therefore the recursion relation for the potential term has the same 
form as ([29|): 

i ^ (30) 

2.4 Justification of the localization approximation; Appli- 
cation to constant fields 

Let us summarize what we have achieved so far. Given the functionals and 
W we can compute the fluctuation propagator via equation (|2^) and the 
derivative of the background field via equation (|2^). This enables us to 
calculate VF*"*"^ and F*"*"^ by means of recursion relations (26) and (|^). At 
least this can be done in principle. Because of the functional nature of these 
equations the calculation is too complicated for numerical purposes. We will 
arrive at a significant simplification by arguing that we may focus attention on 
constant fields. 

Fields $ with very large derivatives have very small probability because of 
the kinetic term in the action. We may therefore assume that $ does not have 
such violent behavior. 

\l/*[<i>](z) is an interpolation of the field <i> from the coarser lattice A'+i to 
the lattice A' . The interpolation is determined by a minimality condition on the 
action. Because of the kinetic term in the action, we expect that the result of 
the interpolation is a smooth function on A*. This means that we may consider 
\l'*[<i>] as nearly constant over distances a^, assuming that the scaling factor 
s — ai+i/tti is big enough 

By locality S'\ztu[5'*[$]] is very nearly zero if z is not within about one lattice 
spacing of w, and it depends only on the field within a neighbourhood of about 
one lattice spacing ai of z and w. Therefore we may approximate 14^* [^E"*] by a 
hermitian operator which involves only the value ^ = 5** (2) of the field at 
one point (see also (|l2|)), 

W'[^']{z,w) « ]^(^\z,w\^\z)) + W\z,w\^'{w))^, (31) 

This justifies the localization approximation for VF*+^. VF* has locality proper- 
ties similar to a Laplacian. 

To determine PF*+^[<i>] we need F*[(/)] for — ^'*[$]. F*(a;, y) depends on the 
field on a domain of diameter about a^+i. Unless $ is constant, the field 4> 
need not be nearly constant over such a big domain. But r*[0] is determined 
by S^"[(j)] = — VF* and we can use the localization approximation for this latter 
quantity in order to determine F* for general field 4> if we need it. 

Looking now at the recursion relation for PF'+^(a;,?/) we see that three dif- 
ferent quantities are involved: 

• The fiuctuation propagator has an exponential decay of one block lattice 
spacing a^+i (see [|l| and figure^). 
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z 



w 



Figure 4: The decay of the fluctuation propagator. 



• The operator VK*(z, and its derivatives are very local by assumption 

and have therefore a range of one lattice spacing Oj. 

• The derivative of the background field \E'*' is also determined by these 
quantities and has exponential decay with decay length aj+i. 

Combining these arguments we see that has similar locality properties 

on the scale of the new lattice spacing Oj+i as had on the scale aj. 

Therefore we will be entitled to make similar approximations on the next 
level as before. 

Now we can perform functional derivatives of for almost constant fields: 

+ W''{ZI,Z2\<I){Z2))S{Z-Z2)) (32) 
W\,^mZuZ2) = ^(W"'{ZI,Z2\^{Z1))5{Z - zi)S{w - Zl) 

+ W"'{zi,Z2mZ2mz - Z2)S{W - Z2)) , 

(33) 

where ' denotes the derivative of W^{zi, -22 10 with respect to the single variable 
^. These formulas arc valid when smeared with test functions in variables z and 
w which are nearly constant over distances aj. 
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Setting the block-spin field constant the fluctuation propagator becomes a 



function of the field value. By virtue of (27) the derivative of the background 
field can be written as a simple derivative: 

*\,[ $](z) = *^'(2, x\^) ^ f(s{z, zi) + / r (z, w\'^)W'{w, zi\^))c^*{zi,x) . 

i i (34) 

For constant block-spin fields, inserting (32) and (33) into the recursion relation 
(26) for W gives therefore a recursion relation for W^: 



\ I 

Zi,Z2&A 

1 1 



(35) 



W'+\x,y\^)= / W\z,w\'i>)'^'"{z,y\^)'^"{w,x\^) 

+ r*(^l,^2|W'(22,;23|$)r(^3,^4|W(24,^l|$) 

{^'^'(Z2, x|¥)*^'(z4, ym + *^'(^2, a;|$)*^'(zi, 2/1$) 
+ *^'(Z3, a;|¥)«'"(^4, y|$) + -^"{zs, x|$)^'^'(zi, y|$)} 

r(zi,Z2|W"(^2,^l|¥) 

{**'(zi, :e|$)*''(zi, y|¥) + *^'(z2, x\^)'f''{z2,y\^)] 

r (Zl, Z2\^)W'\Z2,ZI\^)W"{Z3, Z4|¥)*^'(Z4, x|$) 
{r (Z2, Z3|¥)*^'(Z3, y|¥) r (Z2, Z3|$)**'(Z4, 2/1$) 

+ r (zi, Z3|$)**'(z3, 2/1$) + r(zi, Z3|$)*^'(z4, y|$)}. 

For constant fields we have translation symmetry. Therefore the derivative 
of the potential must be independent of the coordinate: 

=F'(0) = const. (36) 

To obtain the recursion relation of the effective potential we insert the local 
approximation for VF*: 

z 

-i J r(zi,Z2|$)W"(^2,2l|$)(*\.(zi|$)+r,4z2|$)) . (37) 

Zl,Z2 

We simplify this. By assumption of section |l| 

C{x,z) = 1. 
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Differentiating the constraint ( p^ ) 

j C\x',zWAz)^5{x' -x), (38) 

z 

and integrating over x' we see that 

/ ^'\,[$](z) = l. 

This leads to the simpler recursion relation 
F+i'(l) =F'($) 

-i j r(zi,Z2|W(^2,^l|l)(*\x(^l|$)+*\.(^2|^)) . (39) 

One can use this to determine the constraint effective potential, i.e. the 
potential on the last lattice which consists of one point only. 



3 Numerical results 
3.1 Getting started 

We want to calculate the constraint effective potential on the lattice numerically. 
We consider an example. For starting action we take the standard (^^-theory in 
two dimensions, 



S' 



^1 jSi\ 



l0°(z)A(z,u;)0°(«;) + j (i™^0"(.)2 + ^0"(.)4 



2 



(40) 



To calculate W^[^ and we need W^[(i)% W^'[(jP], W°"[if], V^'[(f] and 

r°[<^°] for constant fields (jP — <p. We write (see also (13)) in a symmetric 
form: 

W\z^,z^\<jP) = A(zi,Z2)-(m2 + ^0°(zi)^)<5(zi-Z2) (41) 

= A(zi, Z2) - + ^{4P{zif + <^°(Z2)'))'5(Z1 - Z2X42) 

Therefore we have for W'^ 

W°(zi, Z2IO = A(^i, ^2) - (ml + ye')'5(zi ~ ^2) (43) 

which is not an approximation because equation (|3^) is exact in this case. Now 
it follows for the constant field (p 

W0(zi,^2|^) = A(zi,Z2)-(m2 + ^0')5(zi-z2) (44) 

W°'(zi,Z2|^) = -X4>5{zi-Z2) (45) 
W°"(zi,Z2|^) = -\5{zi-Z2) (46) 
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Similarly, the potential 

= v'A4P] \f^^ = ^0^+ (47) 
3.2 Relations in momentum space 

We have fixed the starting point. Next we want to calculate the recursion 
relations for and . To reduce the amount of work we switch to momen- 
tum space. Because we insert constant fields we have translation symmetry. 
Therefore the Fourier transformation looks quite simple. For the notation see 
appendix 

W'{zi,Z2\4>) = [W\p\^)e'P^''-''^ (48) 



p 

^\4z\^) - J vl/^(p|^)e'f(^--) (49) 
p 

r'(zi,z2|0) = J r*(?,q,/'|0)e*((«+')"i-(9+'')"^) (50) 

l,q,l' 

and similar for the derivatives, which are now simple derivatives with respect 
to 0. F* has only translation symmetry on the block lattice. Therefore it is not 
diagonal in momentum space (see appendix |c|) . 

Transforming the recursion relations into momentum space leads to 

W'+\q\^) = J ^"*iq + l\^)W\q + l\$)-^"{q + l\4>) 
I 

1 



FX^i, g', l2\<l>)T\h,q + g', k\4>W*{q - h + kl^l^Wiq - h + k\^) 

q',ll,... ,li 

^''{q + q' + h\4>) + W^'iq' + h\^)}^'\q' + h\4>) + W'' {q + q' + 1^)} 
+ i / T\h,q'M\^W'*{q-h + h\4>W'{q-h + h\4>) 

{W'(g' + lS) + W'"{<l' + h\^)] (51) 
V\h,q', l2\4>)T\-h + h, 0, Z4l^)^^'*((? + hl^Wiq + I3 + IM 



4 

q',h,h,h 



1 

8 

q',ll,... ,li 



[W''{q' + h\<t>) + W''{q' + l2\<t>)}{w''*{q + ^sl^) + W {h\<t>)] 

and 



r+i'(</.) = r'(,^) - i / r(/,g,ri(A)r'(-/ + /'| 



4 

q,l,l' 



[w^'{q + m + W^'{q + l'\^)]. (52) 
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We need also the formula ( p7[ ) in momentum space: 

+ l\4>) = C"{q + l)+ f r(/, q, I'lm^q + I'^C" (q + I') 

i (53) 

3.3 Algorithm 

We are now able to calculate all quantities needed. The following algorithm 
gives an overview of the procedure: 

0. i = : Given the starting action fix the values of W^{p\(l)j) and y°'(0j) 
for a number of constant field values 4>i, ■ ■ ■ ,4>n ^-^d for all momenta p. 

1. Calculate W'''{p\(l)j) and W^"{p\(j>j) by numerical differentiation. 

2. Determine r'(Z, q, l'\(j>j) with the help of the explicit formula (p^). 

3. Calculate -^"{q + ll^j) by (|3|). 

4. Insert all into recursion relations for iy+^(g|(/)j) and V'^'^^' {(pj). 

5. Increment i and goto 1 until i = N. 

We want to see numerically how this procedure works to produce the correct 
- derivative of the - effective potential. Therefore we wrote three programs to 
compare various methods. A Monte-Carlo program gives us reference values. 
A combined heatbath and metropolis algorithm is used. The heatbath is used 
with the kinetic term of the action to produce a Gaussian distributed random 
value for the fluctuation field value. The following metropolis decides on the 
basis of the whole action if this value is accepted. 

As an alternative we calculated the effective potential perturbatively, using 
the Gawedzki-Kupianen formalism The action S''[(j)] is split into a kinetic 
term i0(w*)~^0 and a potential is the massless free propagator. Again 

the notion of fluctuation fleld C with C*C* =0 is used and we have a background 
field «-*: 

^ = + (^ (54) 

The background field is an interpolation of the block-spin field: 



The interpolation operator A'' is defined as 

X' =wT"(v*+i)-i (55) 
here is the free propagator for the block-spin field: 

^ Cv'C* (56) 
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/ \ 



X >"< >-< X X 

>K >^ X- 



Figure 5: Perturbative expansion of the effective potential. External lines rep- 
resent and internal lines F*. 2-point vertices have the weight mf, 4-point 
vertices have the weight and 6-point vertices have the weight 7^. 



Inserting the split of the field (j) (|5^) into we get for the fluctuation integral 



•^■«c (57) 
where the fluctuation propagator is defined a^ 

The fluctuation integral can be rewritten as 

g-y'+M*] ^g^lfrr'^g-y-m. (58) 

This can be expanded. All graphs in one loop order and up to three fluctuation 
propagators are taken into account. This leads to the graphs in figure ^. 

The scale factor (block length) is the same as for the saddle point approxi- 
mation and we used a polynomial parameterization of the potentials: 

^Inserting (^sj) and one sees the same formula for as in (^). For free fields and 
5"' ^ are identical. 
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-1 ± 




-2 I 1 1 1 1 1 1 

2 2.5 3 3.5 4 4.5 5 

Figure 6: The derivative of the constraint effective potential is plotted against 
the block-spin <f>. The bottom line is the original potential. The next upper line 
is the saddle point approximation and the top line the perturbative calculation. 
The open squares represent the Monte Carlo results. This was calculated on a 
16 X 16 lattice with blocking factor 2, TOq = —0.8 and Aq — 0.4. Note that the 
minimum of the potential is where V' ~ (see figure pO| .) 



3.4 Result 

As a result we get the following picture: 

• For small field values the perturbative calculation is comparable in accu- 
racy to our method. 

• For larger field values the perturbative calculation goes rapidly wrong, 
while the saddle point approximation remains accurate and has the correct 
asymptotics (see figure |6 ) . 

• For the values of coupling constants which we consider, larger blocks are 
better than small blocks (see figure 0). This confirms our expectation from 
the discussion of the localization approximation. 

In figure || the flow of the potential is shown. 
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1.6 1.8 2 2.2 2.4 

* 



Figure 7: The derivative of the constraint effective potential is plotted against 
the block-spin The bottom line is the original potential. The next upper 
line is the saddle point approximation with blocking factor 2. The top line is 
the saddle point approximation with blocking factor 4. It is better than with 
blocking factor 2. The open squares represent the Monte Carlo results. This 
was calculated on a 16 x 16 for tuq = —0.4 and Aq = 0.4. 
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3.3 



3.35 



3.4 
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3.55 



3.6 



Figure 8: The flow of the derivative of the constraint effective potential is plotted 
against the block-spin <I> in the vicinity of the minimum of the potential. The 
lines represent V^' for i — 0,1,2,3,4 from the bottom and the open squares 
the Monte Carlo results. This was calculated on a 16 x 16 lattice with blocking 
factor 2, m-Q = —0.4 and Ao = 0.4. The extension of the lattice is getting smaller 
after every blocking step by factor 2. On the other side the lattice spacing grows 
by factor 2, so that the lattice volume stays constant. No rescaling to the unit 
lattice is performed. 
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4 Improved saddle point method with normal 
ordering 

4.1 Advanced method with normal ordering 

The localization approximation becomes more and more accurate the larger the 
scaling factor (ratio of lattice spacings) s — ai+i/ai. On the other hand, the 
saddle point approximation becomes exact in the limit s i-^- 1 (when the nature 
of the cutoff permits such a limit). It becomes less accurate with increasing 
scaling factor because the phase space for high frequency modes increases. This 
is one reason why we choose to consider a fairly sophisticated improvement of 
the saddle point method. 

The self consistent improvement of the saddle point integration consists of 
two steps. In the first step, the old action S^[(j)] = S^['^^ + C] is self-consistently 
normal ordered in the fluctuation field ; at the same time the background field 
v[/« = \I'*[$] is determined. Symbolically 

S'[(l>]=:T'[<f\C]: (59) 

The precise definition and properties of the normal ordered form T' of the 
action will be described presently. In the second step, the fluctuation field 
integration is performed, neglecting normal products of higher order than second 
in the fluctuation field. The integral is again Gaussian, and the result reads 

The trace is again over the space Ti^ of functions C on A* with vanishing block 
average, C'C* — 0. 

The method of self-consistent normal ordering is very old - see |Q ; some- 
times it is called the Feynman Bogoliubov method. Therefore we will refer to 
our method here as Feynman Bogoliubov approximation. The definition and 
properties of the normal ordered action T* are as follows. 

Let dfir be the normalized Gaussian measure (free field measure) with covari- 
ance (propagator) F. Given the fluctuation propagator F*, the normal ordered 
amplitude T* is defined as 

T'[^\C] = J d^l^'m{OS[^' + C+C] (61) 

= expQAr[*']1^5[vi/» + C+e]|^^0 (62) 

The second formula is valid when S*' is differentiable. The first formula is 
more general, it can be used to define a normal product expansion also for 
discontinuous actions (cp. the section ^. 

We use the following notation for functionals with two arguments: 

— - — C] = A J** C] Al'i'' C] = A.J^'* C] 
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The notation is a little sloppy. Because C^C = 0, the second expression only 
makes sense after smearing with a test function / which obeys C'/ — 0, viz. 



T = 0- 



The background field = 5'*[<i>] will be determined by the condition that the 
expansion of the action in powers of the fiuctuation field has no linear term: 

T'[i'\C]^T[-i''\0] + ^CT',,[^\0]C + ■■■ forCC^O- (63) 
That is the saddle point condition 

/ Ciz)T\4'f'] = if = (64) 

J z 

which is equivalent to 

T',[-f'] = C"X\ (65) 

C¥ = $ . (66) 

A'[<i>](a::) are Lagrange multipliers, x G A'+^. 

Now we specify the fluctuation propagator F* as pseudo inverse of T^-- so 
that 

r[w']T\,p']r[^'] = r[i'*]. (67) 

CT'[-^']^r[¥]C" = 0. (68) 

As in the saddle point approximation the fluctuation propagator satisfies the 
constraints (^), is positive definite on 7i^, and equal to the inverse of T^--, on 



But T';; depends itself on the propagator F'. Therefore (67) leads to a gap 
equation for F'. 

One sees that the background field is the minimum of the normal ordered 
form of the action 5"* on the surface with prescribed block-spin $ = CVj and 
F' is the inverse of the Hessian TV; at this minimum; it depends on ^f'. (The 
Hessian is a quadratic form on the tangent space to the surface.) 

The determination of the fiuctuation field propagator F' involves solving a 
gap equation. The gap equation is obtained by differentiating eq. ( |62| ) twice at 

C = 0. 

T\, = e(i*^'[*l*)5"'[*^+^]|^^Q (69) 
F*[^'*] = T\,p\0]-^ onHc- (70) 

The gap equation is shown in graphical form in figure Figure ^ explains 
the graphical notation. Further derivatives can be applied and we obtain a 
generalization of eq. (p3) as shown in figure nu. 
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Figure 9: Graphical notation for normal ordered amplitudes 



Figure 10: Definition of the normal ordered amplitudes. Further legs generated 
by ^'-derivatives will be marked by primes 




Figure 11: Gap equation for the fluctuation propagator in Feynman- 
Bogoliubov approximation 
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This completes the discussion of the improved saddle point approximation. 
Eq.(^0|) is still a formidable functional recursion formula. We need to make it 
practical by our second approximation of localization. The localization approxi- 
mation has to be made for the second derivative T*;; and not for S*'". In addition 
to the formula for the second derivative 5**" there is a formula which expresses 
the first derivative 5*' of the ordinary action in terms of the first derivative 
T*; of the previous normal ordered action; this furnishes the potential up to an 
additive constant. 

We have a similar formula for the derivative of the background field: 

^\,{z) = (i-rB*r..)c**(^,x). (71) 

where Tb^ is defined as inverse of TV ^ on 7i^. In contrast to it is defined as 
inverse of the mixed derivative T*; ^ which involves derivatives with respect to 
the fluctuation field and the background field. 

Remark: It can be shown with some effort that the two propagators are ac- 
tually equal, i.e. rB[\E'*] = r' [^'] when is the background field. See appendix 

i 



4.2 Gaussian integration 

We wish to evaluate the integral ([^) for S^^^, using the quadratic approximation 
(|9|j63l) for S\ viz. 

S'*[** + C] =r[*\0] + i : CT\.[-i'\0]C ■■ 
We undo the normal ordering again. Dropping arguments ^ = we obtain 

s'[-9' + c] = T'[-9'] + ^CT\.[¥]C - itrr;j**]r (72) 



The trace is over 7i^. By arguments similar to appendix |A.l| the fluctuation 
propagator has the explicit form 

r = T\r^ - T\r^c"{CT\~^c*)-^CT\r^ . (73) 

As in the saddle point approximation there is an alternative formula for the 
fluctuation propagator 

r[*^[$]]= lim {T\.[¥[<^>]] + kC"C')~\ (74) 



We argue that the last term in ([72|) is a constant independent of ^P*. Multiplying 
T*;; from the right of (O) gives 



rT'.. = l-T\. ^C"{CT\. V")~^CV (75) 
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We take the trace of ( |75[ ) 

tr(rrV) til -tr(T\.~' C"{CT\.~'C")-'C 
= tr 540) -tr 6^(0) 
= ^ 1 — ^ 1 = const. q.e.d. 

Z X 

Now we can again evaluate the functional integral (0) in Gaussian approximation 
with the result 

= r[**[$],0] - itr In r (76) 

4.3 Relations between derivatives of the normal ordered 
action 

We derive a relation between the and ^-derivative by using the following well 
known change of covariance lemma for Gaussian measures ||^ . If the covariance 
depends parametrically on a variable s, then 

^/d,r.(0)A[0].l/d,r.(0)^r.A^M, (77) 

where = dTg/ds. 

We denote ^''-derivatives by commas "," and ^-derivatives by semicolons ";" 
as before. Applying the change of covariance lemma to the definition (^ij) of 
the normal ordered amplitude we obtain the relation 

T\^[^'] = T^;,[**] + itrr,^[5''']r.. [«-']. (78) 

After another ^-derivative we get 

T',. = + itrr,,[*']r;.;^[*']. (79) 

4.4 Recursion relation for the first and second derivatives 
of the action 

We differentiate the recursion relation (|7^) to obtain a formula for the first 
derivative of the action 5'*"'""'^. The derivative of the trln is computed in the 
same way as before. The result reads 



(80) 



We re-express the ^''-derivative in the first term using eq. ([Tq). A cancela- 
tion occurs and we obtain the final form of the recursion relation for the first 
derivative of the action. 

S'+\, = J^^\.A^)T\^,[^']. (81) 
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Next we turn to the second derivative. We differentiate eq. ( |8l] ) once more. 
Since T';^ is the ^-derivative of the normal ordered action we can use eq. ( [78| ) 
again to evaluate the derivative of T'^-z with the help of 

T\z[^''] = -r[^']T\. ,^[^'*]r (82) 

Differentiating the factor ^>^^x{z) produces a term 



'f\.y{z)T-.. (83) 
Differentiating the constraint ( |66| ) twice shows that ^xy G T^c because 

C{x',z)-^\Az)^S{x' ~x), (84) 



implies 

' C\x\ z)-^\xyiz) ^ . (85) 



The saddle point condition ( |64| ) implies therefore that expression ( |83| ) vanishes. 
As a result we obtain the recursion formula for the second derivative of the 
action 

S'+\.y = I '^\.{z)T\^,^^^\y{w) (86) 

J Z.W 

*^x(^) (^rV^ - i tr(rr.. ,XT\. ,^)) vl.^_^(^) (87) 

We neglected to indicate arguments ^E"* of propagators and normal ordered ac- 
tion. We see that there is no 1-particle reducible contribution. 

A subtle feature of this formula is the appearance of two different vertices 
in the last term; one involves a ^'^-derivative and the other a ^-derivative to 
create the external leg. This feature is familiar from Schwinger Dyson equations, 
cp. section || below. Neglecting the difference would only neglect a two loop 
correction, though. 

The recursion relations in Feynman Bogoliubov approximation are shown 
in figures ^2] and They involve propagators which are to be determined as 
solutions of the gap equation (^) , shown graphically in figure |ll|. 



4.5 Derivative of the background field 

The recursion relations involve the derivative of the background field. Dif- 
ferentiating the saddle point condition (|6j) we obtain 

( c("')*^=.(^)^^^,.[*'] = (88) 

J z.,w 
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Figure 12: Recursion relation for 5*'+^' in the Feynman Bogoliubov approxima- 
tion 




Figure 13: Recursion relations for the S^^^" in the Feynman Bogoliubov ap- 
proximation. The second formula is obtained from the first with the help of 
eq. (|7|). 
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We claim that this equation and the constraint (84) are solved by 



^\^{z)^C^*{z,x)- rB\z,w)r,^,J'^*{u,x), (89) 

where C** = C"(C*C")^^ as before and Tb* is defined as the inverse of TV . on 
Ti.},. The proof proceeds as in the saddle point section. 

The remark at the end of the last subsection leads to a simplification of ( p9| ) 
because F* can be substituted for Fs'. See appendix^. 

4.6 Discussion of the localization approximation 

There are new aspects of the localization approximation which need to be dis- 
cussed. The consideration of this section assumed that the starting amplitude 
iS** is accurately known. Because of this, it is appropriate to regard the recur- 
sion relations as recursion relations for the normal ordered amplitudes T*, and 
to make the localization approximation for these amplitudes, not for 5*. The lo- 
calization approximation for is justified by appealing to the same arguments 
as before: In computing 5**+^ we only need T'^; for smooth fields. 

Thus we may regard one step in the recursion relations as composed of three 
parts 

1. Localization approximation on for smooth fields . 

2. Computation of 5*+^ from for smooth fields 

3. computation of T*+^ from 5*'+^ 

When proceeding in this way, the result for S*'^^ must be substituted into the 
(gap) equations which determine the normal ordering. In general there result 
graphs which contain propagators from two successive length scales i and i -I- 1 . 
This complication does not arise when only 1-loop graphs are retained. 

There is one further technical complication. In contrast with S^" , T^-- is 
not necessarily exactly invariant under arbitrary translations by lattice vectors 
in A* but only under block lattice translations. Therefore it it is not easy to 
calculate its pseudoinverse (the fiuctuation field propagator) exactly. One can 
proceed as follows, however. 

Because T^-- and F* are determined by a gap equation one can start with 
some guess for T';;, e.g. with the pseudoinverse of 5*" which was used in the 
saddle point approximation. Then one iterates the gap equation. 

The gap equations are here considered for reasons of accuracy and not be- 
cause there are infrared problems. Quite on the contrary, the propagators must 
have exponential decay. Therefore the iteration can be expected to converge 
very fast. 
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5 Schwinger Dyson Equations 



In principle one can improve on the approximations for the fluctuation field 
integral if the action is known. In any case, more accurate methods are 
applicable in the first step, since the starting action S'^ is known by assumption. 

If S' is only known approximately, there can be problems because the cor- 
rection terms might depend on details of the action which are neglected when a 
localization approximation is made. This is the crucial point to discuss. It will 
turn out that inclusion of two loop corrections requires that also the localiza- 
tion approximation is improved to the next order. One needs to consider third 
derivatives of the normal ordered action evaluated at nearly constant fields. 

Here we discuss the evaluation of the fluctuation field integral by solution of 
Schwinger Dyson equations. 

Let r be a positive semidefinite operator on Ti.. Then the normalized Gaus- 
sian measure d/ir on Ti is defined. Its characteristic function is 

J dfir{(l>) expiqcj) = exp ^-^^r^^ . (90) 

If is a function of which is integrable with respect to the Gaussian measure 
c?/ir, it can be written in normal ordered form, as we know. 

Am = (91) 
sm = + (92) 

The normal ordering prescription depends on F. 

If r is strictly positive and bounded away from zero, then B[(j}\ is an entire 
analytic function of (j) and can therefore be expanded in a power series. This 
yields the expansion of A in normal products of (f>- It exists whether or not A 
is continuous or differentiable. This is important when one wishes to deal with 
singular cases such as a discrete Gaussian model, where the action depends on 
the integer part of a real lattice field. 

Here we are interested in cases where F has zero eigenvalues. The 1-dimen- 
sional Dirac measure S{x)dx is the prototype of such a measure; it is not a 
pathological case. In this degenerate case F projects on a subspace He of Ti 
and B[<j) + ^] is an entire analytic function of ^ g Ti-Ci assuming the restriction 
Fc of F to Tic is bounded away from zero. 

The normal ordering commutes with shifts of the field and with derivatives 
with respect to (j). 

We need the following general integration by parts formula 

dAir(0)AMe-^[^l - S[rA] J d^r(0)e-^[^+«l | ^ ^ p. (93) 
This generalizes the well known special case [|j 

/ dMr(</>)0(^)e-^[^l = / r{z,w)-^ I d/ir(0)e-^[^+«l L „ 
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and can be readily proven with the help of eq. (^0|), as follows. 

Proof of the general integration by parts formula: We use the 
Taylor formula in compact form 

e''^F(0) = F(0 + a). (95) 

Since an arbitrary functional -B[0] admits a functional Fourier transform, it 
suffices to prove the formula (|9^) for the special case 



We compute 



d/ir('?!>)A[(/)]e^90+|9rg 
dtir{(f>)A[(f>] : e*«'^ : Q.e.d. 



The general integration by parts formula has many uses. Suppose we define 
expectation values of an interacting theory 

(A) = Z-iydMr(0M[0]e-^[*l, (97) 
dAir(0)e-^[^], (98) 



and the generating function G[$] of amputated Green functions of the interacting 
theory 

e'^ra = / dA*r(0)e-^[^+«l. (99) 



We preferred to introduce an unnormalized generating function, so that G[0] = 
In Z. The normalized generating function is 

The integration by parts formula tells us that the expectation values equal 

(^)=e-«["lS[rA]eGKl|^^^. (100) 

The Schwinger Dyson equation for the generating function G is another appli- 
cation of the integration by parts formula. It is obtained by differentiating the 
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definition of e*^ with respect to ^ and converting the resulting downstairs factor 
j^V into a differential operator. Suppose 



Then 



V[<j>] =: W[ 



(101) 



w. 



If r has zero modes, this is only valid for the directional derivative in Tic 
directions, i.e. when multiplied from the left with F. Additive constants in G 
drop out. 



Let 



r 



j d^l^{C)e-'"^''''^+^^ for C^H, 



(102) 

(103) 
(104) 



We chose to work with an unnormalized generating function for amputated 
Greens functions again; the normalized generating function is 

T^d"^]^^ is the inverse of the restriction of r'[5'] to its range Hq. Dropping an 
additive constant, the exact effective action is 



-tr lnr[*], 



(105) 
(106) 



The second term compensates for the transition from the unnormalized Gaussian 
measure in definition (^) to a normalized one. It involves a trace over 7i^, cp. 

section ||. 

The Schwinger Dyson equation ( 101 ) carries over literally; V, W, G all depend 
parametrically on ^E*. The single power of C which arises from differentiating the 
second term in eq. (102) converts into a differential operator. As a result one 
obtains 



c(^) 

Ti 



^ g-G*[*,«] 



/ r(z,HT^™[vI/,c + ^'i:]e«•[*^«+^l 

Jw oQ 



saw) 



(107) 
o(108) 
(109) 
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Figure 14: Condition on one- and two-point vertices in the Schwinger Dyson 
approach 



The above Schwinger Dyson equations are true for any choice of the back- 
ground field ^' and for any choice of the propagator F'. A great simphfication 
results if they are chosen judiciously. 

It is customary to call the amputated Green functions vertices, n-point 
vertices are obtained as n-th ^-derivatives of ivT* at ^ = 0. They depend para- 
metrically on the background field. In this terminology the first and second 
derivative if'; and K'^-- of with respect to ^ are the one- and two point 
vertices. Only the directional derivatives in directions in Hq are defined. We 
require 

Clif^[^',e = o] - 0, (110) 

Cli^\;[^,C = 0]C2 = 0. (HI) 

at = 5** [$] and for arbitrary Ci : C2 G '^c ■ These conditions are shown graph- 
ically in figure |l^. In the figures, the vertices obtained by differentiating 
are cross hatched. The Schwinger Dyson equations for the n-point functions are 
obtained by differentiating eq. ( |108| ) with respect to ^ at ^ = 0. The resulting 
equations for one, two and three point functions are shown in graphical form in 
figure |5[ 

In principle the vertices can be computed from the Schwinger Dyson equa- 



tions by iteration. As a consequence of conditions ( |110| , |lll[) the equations for 
the one- and two point functions substitute for the earlier equations for the 
background field 5* and fluctuation propagator F\ In zeroth approximation, 
the vertices with three and more legs are equal to the corresponding normal 
ordered amplitudes ("bare vertices"). The iteration proceeds by inserting the 
previous approximation for the vertices on the right hand side of the equations 
for the three-point functions, and similarly in higher ones. 

Because we have an infrared cutoff a^^, we can expect that the iteration 
converges fast. 



The effective action can be computed fr om eq. (105). We determine its 



i>-derivative by differentiating the definition (103) of G*['I', 0]. Using the change 
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Figure 16; Schwinger Dyson Gap equation in 1-loop approximation 



of covariance lemma ( |77| ) again, we obtain 
S 



y\,[vE',C])e-^'[*^';i (112) 



iAr - 

where F = F*[5']. We restrict attention to ^-derivatives. From the definition of 
we get 



(113) 



/ *,.(^) (: : +i : CF^^F^.F^^C : +^trF, 



This is inserted into eq. (112). We convert again the downstairs factors of C 
into differential operators, using the i ntegr ation by parts formula. One term 
cancels against the second term in eq. ( 105 ) and we obtain the final result 

^ - ^ ^A^)Q"-P,o]{z) (114) 



r:.[vI/,e + F^]e^'[*^«+«l 



C = o 



(115) 



The result is shown in graphical notation in figure |r^. The derivative of the 
effective action can be computed when the X*-vertices have been determined 
from the Schwinger Dyson equation. 

We consider next the equation for the background field. Comparing the 
above definition of Q*. with the Schwinger Dyson equation (eq. (108) at ^ = 0) 
for the 1-point function we see that the equation for the background field takes 
the form 



CQ\i^] = for aU C e Wj; 
The Schwinger Dyson equation for the 2-point function implies 

rQ\.r = r. 



(116) 



(117) 



One deduces (in the by now familiar way) the formula for the derivative of the 
background field. 



^\^{z)^C^*{z,x) 



Tb'{z,w)Q\,,u{w)C'*{u,x). 



(118) 
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Figure 17: Recursion relation for S^' in the Schwinger Dyson approach. 



r^* is defined as inverse of the mixed second derivative . 
have 



on TLq. Also we 



(119) 



One may finally differentiate (114) once more to get an equation for the 
second derivative of the effective action. The term involving the second deriva- 
tive of the background field vanishes again because the constraints imply that 

XV e T-L}' (see (Is^). As a resuh 



5<^(x)5<^{y) 



'5C 



(u) 



One can insert ( [t^ ) and eq. (119) for the first term in order to eliminate a 
reference to the first ^—derivative of T*. 

The equations for the first and second derivative of the effective action are 
shown in Figures 17 and 18, Again, there are no 1-particle reducible contribu- 
tions. 

From these equations one can see what are the leading corrections to the 
approximations which were considered in the previous section. To make the 
comparison we take the zeroth loop order of the Schwinger Dyson equations. 
After inserting b) into c) in figure |l5| in zeroth order we get the simplified vertices 
shown in figure |l^. Using these in the recursion relation for the derivatives of 
the effective action and keeping only 1-loop graphs leads to the familiar result 
of the previous section (figures O and n3) . 



5.1 Two loop corrections 

When we make appropriate localization approximations, the Schwinger Dyson 
equations should become integro-differential equations for matrix functions of 
a single variable. We see, however, that the equations involve normal ordered 
amplitudes with more than two high frequency legs. 
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Figure 18: Recursion relation for S""' in the Schwinger Dyson approach. The 
prime denotes a derivation with respect to the background field. 




Figure 19: Schwinger Dyson equations of the first three vertices in zeroth order 
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The Schwinger Dyson equations for the vertices and the formulae for the 
derivatives of the effective action require a separate discussion. 

We see that the two loop terms in the recursion relation for the second 
derivative of the effective action involves normal ordered amplitudes with three 
high-frequency ("hard") legs. They are not known if we only know the second 
derivative of the normal ordered action evaluated at nearly constant field. 

Therefore we need consider third derivatives of normal ordered actions eval- 
uated at nearly constant fields. This introduces quantities {z, w, u\r]) . They 
have exponential decay with the tree distance of z, w, u with decay length not 
bigger than one lattice spacing. 

We will also need recursion relations for these 3-point quantities. They can 
be obtained from eq. (120) by differentiating once more. We will not write the 
result explicitly, the graphs involved are similar to those in the Schwinger Dyson 
equation (figure |l^) for the 3-point vertices, except that the external legs are soft 
ones. It appears that this equation involves a normal ordered amplitude with 
/oitr hard legs (in the last graph), and so we seem to be in trouble again. But this 
is actually not really so. Because the external legs are soft, there are actually 
only two independent hard relative momenta in this graph as in the remaining 
graphs. The arguments of two hard legs which join to different soft external 
legs may therefore be identified, and the knowledge of up to third derivatives of 
the normal ordered action suffices to determine this vertex accurately enough. 

The strategy for the Schwinger Dyson equations is different. There are 
vertices with still more hard legs. We regard the determination of the fiuctuation 
propagator and of the (cross hatched) vertices as one task, whose input is the 
effective action 5'. We propose to make localization approximations for the 
normal ordered actions T% but not for S^. Therefore S*' may be regarded 
as accurately known for arbitrary fields from the previous application of the 
recursion formula. Therefore one may take derivatives of any order to create 
hard legs. One needs to insert the resulting formulae into the Schwinger Dyson 
equations. As a result, one has to evaluate graphs involving propagators at 
two different length scales. This complication arose already in the Feynman 
Bogoliubov approximation in section ^. 

The practical evaluation of all the two loop graphs requires a serious effort 
in high performance computing, unless one is willing to approximate the rather 
complicated exact formulae for the bare propagators and vertices. We are not 
prepared to elaborate on this. In this paper we only discuss the actual evaluation 
of 1-loop graphs. 

6 Relation to other approaches 

6.1 Flow equation for one particle irreducible average ac- 
tions 

In this paper we are chiefiy interested in blocking to discrete "fiows" from lattice 
to lattice because we wish to investigate the accuracy of variants of the method 
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by comparison with Monte Carlo simulations. But the method itself is not 
limited to this situation. We could consider actions which depend on a 
continuously variable cutoff K in place of the discrete length scales Oj, and we 
may leave k finite, or take the limit k i— > cxd at the end. The flow will depend 
on a substitute for the projector 

Rk = kC"C\ (121) 
Differentiating the recursion relation in saddle point approximation, eq. (p[) with 



respect to t = In if, we get from the formula (25) for the fluctuation propagator 



^^S'' = \tr{S^" + Rk)-'^Rk. (122) 

This is Wetterich's flow equation for average actions |Q. The initial condition is 
also the same. We conclude that Wetterich's effective actions can be interpreted 
as Wilson effective actions, considered as functions of the background field. 

This is a remarkable result because Wetterich's average actions are not by 
definition equal to a Wilson effective action. They are one-particle irreducible 
while a computation of a Wilsonian effective action in perturbation theory con- 
tains one particle reducible pieces. A one particle reducible piece occurs also in 
Polchinski's version of a flow equation [ pO[ . 

However, there are no one-particle reducible pieces in the recursion relations 
when the Wilson effective action is considered as a function of the background 
field. For given values of the block-spin, the background field depends not 
only on the block-spin but also on the action. When the background field is 
determined accurately enough (e.g. in Feynman Bogoliubov approximation) no 
one particle reducible graph appears. 



6.2 Perfect actions 

Let us also clarify the relation of our approach with the work of Hasenfratz 
and Niedermayer iQ on perfect actions. They also use a background field 5". 
Let (T : A* I— > A'+ be the map which scales the coordinates of every lattice 
point by a factor s = ai+i/ai. This induces a map of fields a* : i— > W 
defined by cr*$(z) = s'<i>((Tz), where / is determined by the dimension of the 
field. When seeking a perfect action in tree approximation, the irlnF* in the 
recursion relation is neglected, and the fixed point condition which defines a 
perfect action in tree approximation reads 

S\^'[^) = S\a*^). (123) 



6.3 Legendre transforms of higher order 

It was pointed out to us by Yu. Pismak that our method invites the ap- 
plication of the formalism of Legendre transform of second order jlSj, and of 
higher order when higher order corrections are considered. This yields expan- 
sions in skeleton graphs similar to those familiar from renormalization theory 



38 



Q with two-point functions and vertices which are determined as sohitions 
of the Schwinger Dyson equations. 

7 Theories with fermions 

In principle the considerations of this paper apply also to theories with fermions, 
if one knows how to block them. The (5-function in the defining recursion relation 
for the actions has only a symbolic meaning in this case, but the split of 
the integration over Fermi fields into integration over high and low frequencies 
is nevertheless possible in the same way as for Bose fields. But because of the 
different formula for Grassmannian Gaussian integrals, the second term in the 
recursion relation (^) changes sign. In the formulae for normal ordered actions 
and for derivatives of actions, this is taken into account by the familiar rule 

a factor (-1) for every closed loop 

It is known how to block from fermions with the right number of flavors in the 
continuum to Kogut Susskind lattice fermions, and on from there. How to do 
this in a way which preserves locality was first discovered by G. Mai and is 
reviewed in [ pTf ; Bietenholz and Wiese found a similar scheme in their studies 
of perfect actions Q . 

In fermionic theories the improvement of the saddle point method by self- 
consistent normal ordering can be important, see the next section. 



8 How locahty can fail 

Things can go badly wrong when locality of the effective action fails due to a 
bad choice of block-spin. It is therefore essential to monitor whether the method 
can be expected to give reliable results by monitoring the range 

U{z-w)w{zMn 

of the interaction for each value of the constant field which enters into the 
computation of the desired quantity. 

Failures of locality can occur for several reasons which require different reme- 
dies. 

One sees from the recursion relation for the derivatives of the action that 
good locality properties can only be expected if the fluctuation propagator 
T''[z,'w\^) decays with distance (z — w) with decay length no more than one 
block lattice spacing . 



T = lim [ -A + + ^X<f){zf + kC"C ] . (125) 
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Figure 20: Potential and its derivative. The middle region is unstable and 
unphysical. It is separated from the stable region by a metastable region. Our 
method cannot compute the constraint effective potential in the unstable region 
because we assumed that translation invariance is not spontaneously broken 

Following Balaban j^], one can estimate the decay properties by the lowest 
eigenvalue of the operator in brackets with A replaced by the Laplacian with 
Neumann boundary conditions on the block boundary; basically this eigenvalue 
should be strictly positive and of order at least one in units where a^+i = 1. We 
are interested chiefly in (nearly) constant </>; in this case one can also find the 
decay by Fourier analysis cp. ref. The KC"C*-term effectively eliminates 
the zero mode of A^ (constants) and the next eigenvalue of A^ has the desired 
magnitude. Therefore things go wrong when m^ + ^A0(2;)^ becomes too strongly 
negative. This can happen when « and is too strongly negative - i.e. 
near the top of a Mexican hat which is too high, cp. figure 

When is negative enough, something still more terrible happens. Nu- 
merical work revealed that the minimal action for constant block-spin zero is 
reached for background fields which have periodic domain walls of alternating 
slope. In other words, the auxiliary theory of section |^ shows spontaneous 
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breaking of translation symmetry by block lattice vectors. 

At intermediate RG-steps one could try to remedy this by choosing block- 
spins of fixed length - i.e. integrating out the modulus of the field - to obtain 
nonlinear cr-model type effective actions. 

But anyway we cannot calculate the dependence of the constraint effective 
potential on the magnetization in the unstable region (figure |2^) because the 
assumption of no spontaneous breaking of translational symmetry is physically 
wrong there. 

Nonlocalities can appear in another more interesting way which is relevant 
for studies of a dynamical Higgs mechanism, in particular superconductivity. 
The remedy in this case is the introduction of a composite Bose field as a block- 
spin. In superconductivity it represents Cooper pairs. This mechanism was 
studied in some detail by Grabowski . 

Nonlocalities of 3-point functions will also induce nonlocalities of S''". The 
formal solution of the Schwinger Dyson equation for the 3-point function in- 
volves (among others) the famous chain-of-bubble diagrams. The fluctuation 
propagators in these diagrams have exponential decay. Nevertheless the sum of 
these diagrams can fail to have the desired locality properties due to a pole in 
momentum space below the UV-cutoff a~|_\ of the new action. This occurs in 
fermionic theories when the coupling gets strong enough. In the 2-dimensional 
Gross Neveu model and in models of superconductivity, one is always driven 
into this domain by the renormalization group flow. When this happens, one is 
forced to introduce a composite block-spin ||ll| . 

9 Motivation for the improvement of the saddle 
point method 

In a (^^-theory the simple saddle point method works quite well. But this is due 
to the simple form of the kinetic term in this theory. 

Consider instead an action for an n-component field (j) in two dimensions of 
the form 



where / is a coupling constant. In the continuum, (5(0) is quadratically diver- 
gent, on the lattice (5(0) = a^^. 

This action comes from the 0(n + l)-symmetric nonlinear cr-model in stere- 
ographic coordinates; the term proportional (5(0) comes from the measure on 
the sphere. But let us forget where the action came from. 

Suppose we wish to block from the continuum to the lattice. For simplicity, 
consider first what happens if we apply a simple saddle point treatment to the 
whole theory. Expanding the action up to terms quadratic in the field, we get 





(127) 
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Evidently we met with disaster. There is a quadratically divergent coefRcient 
in a quadratic action. The same problem appears if we do the high frequency 
integral to compute an effective lattice action from the continuum action. 

If we normal order first before expanding in the field, the situation is dif- 
ferent. Let us imagine a finite volume and normal order with respect to the 
massless free propagator vcb without its zero mode. Normal ordering the lead- 
ing term /(V0)^(/)^ produces a quadratic term ~ f AvcbiO)4>'^ = /(5(0)(/)^ which 
cancels against the quadratically divergent term from the measure. 

The disease of the simple saddle point method comes from its lack of covari- 
ance under field reparameterization. This can lead to catastrophic violations of 
Ward identities in theories with symmetries. f\ Since we shall want to apply our 
method to such theories in the future, it was essential to go beyond a simple 
saddle point approximation. 

Another advantage of the normal ordered version of our method (section 
^ consists in the fact that it enlarges the class of theories which can be dealt 
with to theories in which the field assumes a discrete set of values, such as the 
discrete Gaussian model, Ising and Potts models etc. This is true because any 
action can be expanded in normal ordered products of the fields. 

10 Models whose fields assume discrete values 

Here we wish to explain how normal ordering can help to deal with models whose 
fields assume a discrete set of values. Only the main idea will be presented; 
detailed numerical studies remain to be done. 

Consider as an example the discrete Gaussian Model in two or three dimen- 
sions. It lives on a lattice A°. The field assumes integer multiples of 2tt as its 
values. The starting action is 



In two dimensions, this model is related by a duality transformation to the plane 
rotator model (with Villain action |^) which has a Kosterlitz Thouless phase 
transition. In three dimensions it is the dual transform of a U(l)-lattice gauge 
theory. The renormalization group flow of this model is well understood and 
lead to a rigorous proof that the three dimensional U(l)-lattice gauge theory 
shows linear confinement for arbitrary values of the coupling constant . The 
overall factor in the action is written as because the duality transformation 
interchanges high and low temperatures. 

The basic idea is to regard the discrete field as a function of a continuous 
field (j). One regards the action as a function of this continuous field. One normal 
orders it in a self-consistent way. This furnishes at the same time a split into 
a free action which is quadratic in and an interaction. In favorable cases, a 
good approximation to the self-consistent split can be guessed a priori. 

■^It is known that Ward identities guarantee a correct cancellation of divergences in the 
nonlinear ct- model |6|. 




(128) 
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The normal ordered action is an entire function of the field and can be 
expanded in powers of the field components which one wishes to integrate out. 

Let N{S^) be the periodic function of ^ S R with period 27r/3~^/^ which is 
defined by 

iV(0 = S for - 7r/3-i/2 < ^ < 7^(3-^/^. 
We note that N{£^) is small for all values of 1^ if /? is large. Now we may write 

n{z) = [^(^) _ 

We may substitute an integration over 4>{z) for the sum over n{z). The action 
takes the form 

S'W = ljN,cj^{zr + V,'icp) + V,'i<P), (129) 
yO(0) - / iV(0(z))A0(z) (130) 



^2°(0) = I I N.Nicj^iz)^. (131) 



Given the propagator T with respect to which we want to normal order, the 
normal ordered form of the action is determined by Gaussian integration ( |6l[ ) . 
In the case at hand, the result can be computed by Fourier transformation. 
(There are also more direct methods which exploit the fact that the Lagrangian 
depends only on the field at two lattice points |]l6| ). 

The periodic function N((p) of </? e R admits a Fourier expansion 

iV(^) ^ ^me''"'3'''^ (132) 

mez\{o} 

Nm = -r'/'(-l)™. (133) 
m 

We wish to normal order in the fluctuation field ( which is to be integrated 
out in one RG-step. We decompose as usual (/) = ^ + 5" = Now we can 

write the three terms in the action in normal ordered form. The quadratic term 
is trivial and remains quadratic. We use partial integration and the formula ( pO| ) 
for the characteristic function of a Gaussian measure to compute the normal 
ordered form of V^. 



dfiriO I (A^{z)Nii^ + 0{z))+ I Ar{w,z)-^N{i<l> + 0{z)) 



^0 (134) 
A0(z)sin(m/3^/20(z)) + /3i/2m / AT{w,z)cos{m[3^/^(t){z)) 

J W 

- «' + C- 
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The second term V2 can be treated in a similar fashion. First the lattice deriva- 
tives must be written as differences of two terms. After that, the Fourier expan- 
sions are inserted and the Gaussian integration can be performed as before. We 
will not give the somewhat complicated result in full, but only an approximation 
to the resulting effective action which is valid when f3 is large. Although it is not 
quantitatively correct when f3 is not large, it will suffice to discuss qualitative 
features for orientation. The V2'-term is small compared to the T^j^-term in this 
limit. 

The crucial aspect of the model is the breaking of the symmetry under field 
translations 

(t>{z) (l){z) + 27r/3"i/2a. 

The free action is symmetric under translations with a G R, but the interaction 
breaks the symmetry down to translations with a G Z. 

When (3 is large, the expressions exp (— /3m^r(z, z)/2) are exponentially 
small for m ^ 0, and the leading symmetry breaking terms are those with 
ni ~ ±1. We neglect to write the others, and also the small symmetry breaking 
corrections to the generalized kinetic term . 

The effective action after one step is given by eq. (60). It depends on 
the normal ordered form of S*" at C = 0, on the fluctuation propagator and on 
the background field as we know. It takes the following form in d dimensions 

S\'^>) = ^ Q[V*(zr + A(z)cos(/?i/2vI/(z))-ilnr(z,z) + ...^(135) 
A(z) = 2e-'^r("'")/2 / Ar(w,z). (136) 



* = (137) 
F = F°[^']. (138) 

For a translation invariant propagator F, A(z) would be constant. The propa- 
gator depends on the dimension d. 

Below we will briefly discuss some qualitative implications on the basis of 
this approximate expression, neglecting a possible ^'-dependence of F. ^ Quan- 
titative calculations could start from the exact normal ordered expression for 
5°. The exact 1-loop effective action after one step is still given by eq. (|60|), 
which can be differentiated. A localization approximation is made after the first 
RG-step, and the calculation proceeds in the usual way from there on. 

The fluctuation propagator depends on an infrared cutoff M — a^^. For 
orientation, we may imagine that the cutoff is lowered by as much as we please 
in one step; we discuss what happens if M 1— > 0. 

In d — 3 dimensions, the propagator F has a limit when M 1-^ 0. Therefore, 
also A has a limit. However, to judge the RG-flow, we need to rescale to the 
unit lattice. Therefore we need to consider the dimensionless quantity afX. It 
goes to inflnity. Therefore the cosine shaped potential wells become inflnitely 

•^The standard free \I/-independent fluctuation propagator of Kupiaincn and Gawodzki ^ 
gives a good a priori guess for tlie self-consistent split in this particular model. 
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high in the infrared hmit, and there will always be spontaneous breaking of the 
Z-symmetry with a nonvanishing surface tension. There is also a curvature at 
the minimum which produces a mass. In the dual picture, the surface tension 
becomes the string tension of the ?7(l)-lattice gauge theory. 

In d = 2 dimensions, F diverges logarithmically when M i— > 0. The asymp- 
totic behavior of afX as M i-^ depends then on /3. There will be a critical 
value /3c- Its precise value depends on details which were neglected in ( |138| ). For 
f3 < Pc, afA tends to infinity, and for /3 > /3c it tends to zero. This corresponds 
to the two phases of the plane rotator which are separated by the Kosterless- 
Thouless phase transition. The Z-symmetry is either spontaneously broken or 
enhanced to R. This phenomenon of symmetry enhancement was discovered 
long ago by Frohlich and Spencer 0. 

It might also be of interest to consider the dual transform of the 4-dimensional 
[/(l)-lattice gauge theory which is a Z-gauge theory. In principle it could be 
treated with our method. Discussions on how to deal with gauge theories are 
found in glig. 
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A Formulas for saddle point approximation 
A.l Form of F 

We are looking for the propagator which corresponds to the quadratic term 
iCS""[*i*]]C in (21). The fluctuation fields C' obey the constraint C'C = 0. 
Leifl 



It obeys the constraint CT' = 0. Therefore F* : W H^. We now want to 
show that 



This means that 

The first equation is equivalent to 



Applying S^"C to (139) from the right confirms this. 

The second equation is true because it is the adjoint of the first. 

*S"' has in general no zero modes. If there are some zero modes (e.g. for </< = 0) then they 
can be ehminated with the help of projection operators before building the inverse S'" 
This is well known (see 
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A. 2 Proof of the relation for the background field 



Differentiating the defining equation ( |1S| ) for the background field by the chain 
rule, we obtain 

S'"[^i]'qi'' ^C"X'' (140) 
for some A*' e Muhiplying ( |l40| ) with S'"~^ 

^ S"''^C"X" (141) 

and C gives 

By virtue of C'*'' = 1 it follows 

1 = {CS'"~^C*)X'' 

Therefore 

{CS'"~^C*)-'^ = X" (142) 
We have now one possible form of ^E**': 

^ 5''"-V"(C*S''"~V")-^ (143) 
Multiplying in ( |13S| ) with S^" from right and using (142) gives 



= 1 - S'" ^C"X''C 



Inserting (141) leads to 

rS'" = 1 - **'C\ (144) 
By multiplying C* from the right we get the result 

= (1 -rS"")C**. (145) 

A. 3 Derivative of 



Multiplying (144) with from the right and using CT' — gives 

-pi gill^i pz 

Differentiation leads to 

pi/ pi/ ^i//pi _|_ pi giin^i pi gill-^il 



S'"T^ and F^S"*" can be replaced now by (|l4j). Using CT'' = leads to the 
result 

pi/ pi 

or expressed with W'' 

F\^ = F*Vl^%r* . (146) 
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A. 4 Determination of the background field 



An equation for was already derived in A^, It involves quantities 

which are obtained together with the propagator. Expressed in terms of it 
looks 

= {i + rw')d*. 

To determine the recursion relation of we need also the second derivative. 
Using the formula for the derivative of the fluctuation propagator (146) we find 



T\z, Zi) W\UZ1,Z2) *\j,[$]H *\x[$](^2). 



A. 5 Recursion relation for 



(147) 



To get W^'+i out of (24) we need the second derivative of S^^^. 
Consider first 

^trlnr[vI/'[$]] ^trrmn s^^^Wm] 



where we have used ( |l4q ) and Vl^',2[5'*] or denote the functional deriva- 

tives 

S S 
We differentiate once more to obtain the second derivative and omit the func- 
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tional argument \E'*[$]: 
52 



■tr In r 



Z,'W.ZiEA'^ 



(148) 



J r (Zl , Z2)W\,.^ {Z2 , [$] {Z)^\y [$] {W) 



Z,'W.Zi,Z2GA^ 

+ J r\zi,Z2)W\,{z2,Zl)^\,y[^]{z) 

Z,Zi,Z2SA* 

Now we consider S'*[^' *[<&]]. By the chain rule 
52 



(5$(x)(5$(y) 



S'*[^''[$]] = ^'\^[$]* S"'[^'] '^\y[^] + S"[^'] '^\a:y[^ 



(149) 



By the saddle point condition, S'"[^''] — C"X\ Therefore the last term looks in 
coordinates 



xiGA»+i 



Inserting (147) shows that this term vanishes because of F'C" — 0. Therefore 



W'+'[mx,y) = 



1 6^ 



2 S<^{x)6<S>{y) 



tr In r 



Using (148) and (147) leads to the result (26). 



B Simplification of mixed derivatives 

In the case of Feynman Bogoliubov approximation we defined the propagator 
r' for arbitrary fields 4> as inverse of S"" on and derived some formulas in a 
general manner. Afterwards we used special fields namely the background field 
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which depends on the block-spin. Now we regard it only as a function of the 
background field. It is convenient to extend this definition to arbitrary fields by 
specifying that it is independent of the high frequency components of the field, 



for C e Ti-c- 

It is important that the formula for the derivative of the propagator 

6 



(150) 

(151) 
(152) 



is then only valid in the restricted sense that it can be used for evaluating 
^-derivatives of the propagator. 

With this specification of the propagator, the definition (62) of the normal 
ordered amplitude implies that 



(153) 



for C £ T^c- other words, T['I'*,^] depends actually only on 'J* + ^ so long 
as ^ e H};. 

We find it convenient, however, to extend the definition (61) of T['i>^,£] to 
arbitrary ^. Eq. (153) remains valid, of course, but Tf^E**,^] is no longer only a 
function of the sum of its arguments. 

When evaluating $ derivatives, eq. (151) can be used to express the deriva- 
tive of the propagator. On the other hand, ^'^-derivatives in directions 
are the same a s ^-d erivatives; the last term does not contribute in this case 
because of eq. (150). 

Also the equation (|9|) for the derivative of the background field simplifies, 
i.e. Tb' = r': 



It fulfills also (|88|) and the constraint (^4|). One needs the identity 



r(*,z) 



(5«''(z) 



(154) 



(155) 



and ^-derivatives of T* agree in Ti^-directions, and the range of is TQ^. 
Therefore we may convert the ^''-derivative into a ^-derivative. The assertion 
follows now from the gap equation (|6^). 



C Momentum space representation 

Let A* be a ZJ-dimenional position space lattice with lattice spacing and length 
Li = Niai. The corresponding dual lattice A* has a lattice spacing di = 27: / Li 
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and length Li = ^-n jai. The block factor between two successive lattices shall be 
denoted by s = ai+i/ai. It is convenient to introduce an additonal momentum 
space lattice A*'*"'"-'^ with spacing a^^i+i — Li and length = sLi. Then any 

p € A' can be decomposed into p — q + I with q e A*+^ and I e A*'*+^. The 
definitions 

imply the split 

(157) 



p q I 



We now want to derive the momentum space representation ( pO| ) of the 
fluctuation propagator. First note that any linear operator L : W ^ W with 
L{z, z') = L{z — z') can be written as 

L{z - z') = / i(p)e'P(^"^'^ (158) 



Examples are W^{z, z'\(j)) in the saddle point approximation and the laplacian 
A*(z, z'). In the latter case one gets 

A\p)^-p' , p, = lsin(^). (159) 

The representation (158) also holds for the averaging operator C^{x, z) since 
x&W+^ CW with 

^J-^ Vssm(p^aj/2) 7 

The fluctuation propagator on the other hand is only invariant under 
translations F'(z, z') = F'(z — x, z' — x) with x £ A*+^ and therefore has the 
more complicated Fourier representation 

F*(z,z')= j r\l,qJ')e'^''+''>'-'^"+^">'' (161) 

l,q,l' 



with 



ril,q,l') = v\q + l)Sij, 

-v\q + l)C"{q + l)u'~\q)Ciq + l'y{q + I') (162) 
u\q) = [ C\q + l)v\q + l)C"{q + l) (163) 



v\p) := W'{^)-\p) (164) 
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